Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 17 August 2010 (MN I^TeX style file v2.2) 



Double-double radio galaxies: further insights into the 
formation of the radio structures 

C. Brocksopp 1 , C.R. Kaiser 2 , A. P. Schoenmakers 3 and A.G. de Bruyn 3 ' 4 

1 Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, Surrey, RH5 6NT, UK 

2 School of Physics & Astronomy, University of Southampton, Southampton SO 17 IB J, UK 

3 Stichting ASTRON, PO Box 2, 7990 AA Dwingeloo, The Netherlands 

4 Kapteyn Institute, PO Box 800, 9700 AV Groningen, The Netherlands 



17 August 2010 



ABSTRACT 

Double-double radio galaxies (DDRGs) offer a unique opportunity for us to study 
multiple episodes of jet activity in large-scale radio sources. We use radio data from 
the Very Large Array and the literature to model two DDRGs, B 1450+333 and 
B1834+620, in terms of their dynamical evolution. We find that the standard Fanaroff- 
Riley II model is able to explain the properties of the two outer lobes of each source, 
whereby the lobes are formed by ram-pressure balance of a shock at the end of the 
jet with the surrounding medium. The inner pairs of lobes, however, are not well- 
described by the standard model. Instead we interpret the inner lobes as arising from 
the emission of relativistic electrons within the outer lobes, which are compressed and 
re-accelerated by the bow-shock in front of the restarted jets and within the outer 
lobes. The predicted rapid progression of the inner lobes through the outer lobes re- 
quires the eventual development of a hot spot at the edge of the outer lobe, causing 
the DDRG ultimately to resemble a standard Fanar off- Riley II radio galaxy. This may 
suggest that DDRGs are a brief, yet normal, phase of the evolution of large-scale radio 
galaxies. 

Key words: ISM: jets and outflows — galaxies: active — galaxies:individual: 
B1450+333, B1834+620 — galaxies: jets 



1 INTRODUCTION 



Fanar off- Riley type II radio galaxies (FRII; | Fanar off &; Ri 
|ley 1974) are well-known for the ejection of highly colli- 
mated jets from a central nucleus. These jets subsequently 
inflate a pair of radio synchrotron-emitting lobes where the 
jets are decelerated by the ram-pressure of the intergalactic 
medium (IGM). "Double-double radio galaxies" (DDRGs) 
form a sub-class of FRII objects in which there are mul- 
tiple pairs of lobes, all strongly coaligned to within a few 
degrees of each other. There are currently ~ 12-17 such sys- 



tems known, identified and studied by [Schoenmakers et al. 
(2000a ), |Schoenmakers et al. ( 2000b|, |Kaiser, Schoenmak- 
ers & Rottgering (2000), Konar et al. (2006), Saikia, Konar, 



Kulkarni|([2006|,|Marecki fe Szablewski| (|2009|, Machalski, 



Jamrozy & Konar (2010) and Saikia Jamrozy| ( 2009; and 
references therein] )) . Recently the first "double-double radio 
quasar", 4C 02.27 has been discovered (Jamrozy, Saikia, & 
Konar|2009| ), showing that the properties of DDRG are not 
purely phenomena of giant radio sources and may even be a 
common, but short-lived, phase of active galaxy evolution. 
The second, inner pair of lobes in each DDRG/Q was 



interpreted as evidence for a second episode of jet activity, 
within the remnant cocoon of the original jet, as opposed 
to knots in an underlying jet. Support for this conclusion 
was increased significantly by the discovery by Brocksopp 
al. (2007) of a third pair of closely aligned lobes in the 
first known "triple-double radio galaxy", B0925+420. Such 
objects therefore give us rare and invaluable opportunities 
to study the duty-cycles of large-scale radio galaxies, despite 
their episodes of jet activity taking place on timescales of 
millions of years. 



Brocksopp et al. ( 2007 ) modelled the three pairs of lobes 



of B0925+420 in terms of their dynamical evolution in order 
to determine the permitted ranges of jet power as a function 
of age of the source and density of the surrounding medium, 
thereby testing the hypothesis that the subsequent pairs of 
lobes were indeed forming within the cocoon of the original 
jet. They found that the formation of the inner pair of lobes 
was indeed consistent with the outer pair having been dis- 
placed buoyantly by the ambient medium. The middle lobes 
were more problematic, requiring higher densities than those 



found within the outer lobes. An alternative model (Clarke 
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Table 1. Angular resolutions, position angles (PA) and la noise levels in each of the new VLA images. 



Frequency 


A-array Res. (PA) 


A-array la 


Time on 


B-array Res. (PA) 


B-array la 


Time on 


(GHz) 


("J ( u ) 


(fiJy beam -1 ) 


Source (mins.) 


n ( u ) 


(fiJy beam -1 ) 


Source (mins.) 


B1450+333 


1.40 


1.7 x 1.5 (-12.4) 


217.0 


236 


6.6x4.9(74.8) 


214.0 


55 


4.86 


0.6 x 0.5 (-21.9) 


28.9 


117 


1.6 x 1.5 (-85.9) 


38.5 


88 


8.46 


0.3x0.2 (-7.3) 


29.1 


117 


0.9x0.8 (-72.8) 


29.9 


89 


B1834+620 


1.40 


1.9 x 1.5 (-8.33) 


212.0 


249 


5.9x4.9(17.09) 


556.0 


62 


4.86 


0.5 x 0.3(13.89) 


44.1 


116 


2.0 x 1.5 (-60.8) 


64.6 


59 


8.46 


0.3x0.2 (-19.4) 


27.6 


98 


1.2x0.9 (-60.8) 


41.1 


57 



& Burns 1991) was suggested, interpreting the inner and 
middle lobes as jet-driven bow-shocks within the outer lobes. 

In order to further fine-tune the mechanisms by which 
the observed characteristics of the lobes can be reproduced, 
it is necessary to model more lobes from other DDRG sys- 
tems. In this paper we use Very Large Array radio data for 
two additional systems, B1450+333 and B1834+620, which 



have been the subjects of detailed studies previously ( |Konar 
et al.|2006||Schoenmakers et al.|2000b] . We summarise pre- 
vious results from the literature in the remainder of Section 
1, we present the VLA observations in Section 2, the result- 
ing total intensity, spectral index and polarisation images 
in Sections 3 and 4 for B 1450+333 and B 1834+620 respec- 
tively and then use the results as inputs for the model in 
Sections 5-7. We note that the modelling sections assume a 
WMAP cos mology, with H = 71kms~ 1 Mpc~ 1 , Q m = 0.27, 
Ov = 0.73 ( |Spergel et al.|2003 ). 



1.1 B1450+333 

B1450+333 (=J1453+3308 ) was identified as a DDRG by 
Schoenmakers et al. (2000a]) using data from Faint Images 
of the Radio Sky at 20-cm (FIRST: |Becker et al.|1995t, the 



NRAO VLA Sky Survey (NVSS: [Condon et al. 1998) and 



the Westerbork Northern Sky Survey (WENSS: jRengelink 
et al.|1997 ). It consists of four radio lobes distributed about 
a central core and only the southern, outer lobe shows ten- 
tative evidence for a weak hotspot. The southern lobes are 
weaker and more extended than the northern lobes and the 
outer lobes display an extended morphology perpendicular 



to, as well as along, the jet axis. Konar et al. ( |2006| studied 
this object further, using a wide range of observing frequen- 
cies and resolutions. They modeled the spectra of the lobes 
and determined spectral ages for the northern and south- 
ern outer lobes of ~ 47 and 58 Myr respectively, implying a 
mean separation velocity of 0.036c; the synchrotron age of 
the inner lobes was ~ 2 Myr, implying a velocity of ~ 0.1c. 



B 1450+333 lies at a redshift of 0.249 ( jSchoenmakers et al. 
2000a). 



1.2 B1834+620 



B 1834+620 w a s similar ly identified as a DDRG by |Schoen-| 
makers et ah] ( |2000a|b ) using data from the WENSS and 
NVSS. It too consists of four radio lobes about a central 
core, with hotspots present at the leading edges of the north- 
ern outer and southern inner lobes. The authors used the 



presence of these hotspots to argue in favour of the inner 
lobes having rapid advance velocities in a low-density envi- 
ronment; such a conclusion was supported by the discrep- 
ancy between the large optical emission line luminosity of 
the host galaxy and the relatively low radio power of the 
inner lobes, as well as by the model of Kaiser et al. (2000). 



B1834+620 lies at a redshift of 0.519 (Schoenmakers et al 



2000b) 



2 OBSERVATIONS 

Radio observations of B 1450+333 and B 1834+620 were ob- 
tained using the VLA on 2000 October 20-21, when the ar- 
ray was in the highest resolution A-configuration, and 2001 
April 2, when the array was in the B-configuration. Data 
were obtained at three frequencies - 1.40-GHz (L-band), 
4.86-GHz (C-band) and 8.46-GHz (X-band) - in each ar- 
ray, as for our previous observing campaign for B0925+420 
(Brocksopp et al. (2007)); IFs at two adjacent frequencies 
were used in each continuum observation. In order to re- 
duce potential bandwidth-smearing in the most extended 
image, the 1.40-GHz B-array observation of B1450+333 was 
obtained in the PA/PB spectral line observing mode. This 
gave us 8 channels of 3.125 MHz each, resulting in a 25- 
MHz bandwidth instead of the more usual, continuum band- 
width totalling 50 MHz. Full polarisation information was 
also recorded. 3C 286 (JVAS J1331+3030) was used as the 
flux and polarisation calibrator while JVAS J1416+3444 
and JVAS J 1849+670 were used for phase calibration of 
B 1450+333 and B 1834+620 respectively and to ensure par- 
allactic angle coverage. 

The 2000/2001 VLA data were reduced using standard 
flagging, calibration, cleaning and imaging routines within 
MPS (version 31DEC04). The flux densities of 3C 286 were 
obtained using the formulae of |Baars et al. (1977) but with 



the revised coefficients of Rick Perley, as is the default option 
in the setjy routine. Additional analysis routines within 
MPS were used to measure lobe parameters and create spec- 
tral index maps. We list the angular resolutions and sensi- 
tivities of the images at all three frequencies and in both 
array-configurations in Table [I] 

We note that these VLA observations of B 1450+333 



have been presented by Konar et al. (2006) and so we do 



not repeat their work here. The B1834+620 observations 
are previously unpublished, although a detailed study of 
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Figure 1. Spectral index map created from the VLA data for 
the inner lobes of B1450+333. We used the B-array images at 
4.86-GHz so as to ensure that no emission was resolved out and 
the A-array image at 1.40-GHz so as to match the resolutions as 
closely as possible. The higher-resolution image was reconvolved 
with the beam of the lower-resolution image to give an angular 
resolution = 1.7 x 1.5". The core appears to have a relatively flat 
spectrum and the lobes a steeper spectrum with a predominantly 
in the range —0.5 to —1.5 



this source has also been presented by Schoenmakers et al. 
(2000b). 



3 RESULTS - B1450+333 

B 1450+333 was detected at each of the six fre- 
quency/resolution combinations. The B-array images clearly 
resolve the core and inner lobes, although the outer lobes are 
only detected in the B-array images at 1.40-GHz and much 
of the extended emission is resolved out in the A-array maps. 
The higher-frequency images hint at an extended morphol- 
ogy to the inner lobes, approximately aligned along the axis 
of the large-scale source. The morphology of this object has 
been studied in further detail by Konar et al. (2006) and 
images at various frequencies can be seen there. 

We used the VLA images to measure the flux, length 
and width of each lobe, typically using the lower-resolution 
images, and to create spectral index images (Fig.[lJ; we de- 
fine the spectral index, a, in terms of the S u oc v a , where 
S u is the flux density at frequency v. The spectral index 
images use B-array data at the two higher frequencies, so 
as to minimise the amount of emission resolved out, and A- 
array data at 1.40-GHz so as to match the resolutions as 
closely as possible. The lobes have a predominently in the 
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Figure 2. Fractional linear polarisation (FP) images of 
B1450+333 using the VLA:B-array data only. FP is represented 
by the greyscale, which ranges from 3a to the map peak in 
each case. The "wedge" indicates the greyscale in the main im- 
age only. The main image is at 1.40-GHz and the inset plots of 
the inner lobes are at 4.86-GHz. The vectors indicate the direc- 
tion of the electric field, following correction for RM (=14 rad 
m — 2 ), and their length indicates the linear polarisation (1" = 
4.81 x 10 -5 Jy beam -1 ). The outer lobes are very strongly po- 
larised - up to ~ 60% in places, particularly around the edge of 
the northern lobe - while the rest of the lobes are still typically 
polarised up to ~ 20%. The vectors show a complicated structure 
to the electric field, particularly in the outer lobes and we discuss 
this further in the text. 



range —0.5 to —1.5. whereas the core is flatter and/or in- 
verted. There appear to be some localised regions of steeper 
spectrum, particularly at higher frequencies. These values 
of a are consistent with those calculated for e.g. B0925+420 
(Brockso pp et al.|2007| } and the lobes of more typical FRII 
radio galaxies, confirming that all the flux has been recov- 
ered. We use these measurements as input parameters for 
the modeling in Section 7. 

Using B-array polarisation images at all three frequen- 
cies, we produced a rotation measure (RM) image. Using 
the CURVAL routine within AIPS, we measured typical pixel 



values of < 10 rad m , consistent with that of Konar et al 



2006, although occasional higher pixel values resulted in a 
mean (using the IMSTAT routine within AIPS) of 14 rad m -2 . 
Using this higher value, we derived corrections to the posi- 
tion angle of polarisation of 32 and 3° for 1.40- and 4.86-GHz 
respectively. Fig. [2] shows the resultant fractional polarisa- 
tion images with electric vectors corrected for RM; the main 
image is at 1.40-GHz so as to show the outer lobe struc- 
ture (and at higher resolution than the polarisation images 
presented by Konar et al.| |2006) and the inset plots of the 
inner lobes are at 4.86-GHz. The electric field vectors are 
roughly perpendicular to the jet direction for most parts 
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Figure 3. VLA images of B1834+620 showing the two pairs of lobes. All contours are at -3, 3, 6, 12, 24, 48, 96, 192xcr. LHS: the 1.40-GHz 
B-array data showing extended emission along the jet axis O = 556 fjjy beam -1 ). Centre: 4.86 GHz B -array data showing discrete lobes 
without extension (a = 64.6 /iJy beam -1 ). Right: 4.86-GHz A-array image showing the inner lobes only (a = 44.1 /iJybeam -1 ). 
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Figure 4. Spectral index maps created for the inner and outer 
lobes of B1834+620; we used the A-array image for 1.4-GHz and 
the B-array image for 4.86-GHz so as to ensure that no emission 
was resolved out at the higher frequencies. The higher-resolution 
image was reconvolved with the beam of the lower-resolution im- 
age, resulting in an angular resolution = 1.9 x 1.5". The core 
appears to have a relatively flat spectrum and the lobes a steeper 
spectrum with a in the range —0.5 to —2. Inset plots show the 
inner lobes at higher magnification. 



of the inner lobes, as is typical for lobes in standard radio 
galaxies. However there are also significant disruptions to 
an ordered field, particularly in the centre of the northern 
inner lobe. The magnetic field becomes perpendicular to the 
source axis in the outer regions of (in particular) the north- 
ern outer lobe. This may be caused by compression due to 
lobe expansion, as also suggested for the middle lobes of 
B0925+420 flBrocksopp et al.|[2007t , which may seem sur- 
prising for a lobe which is no longer supplied with additional 
energy (see Section 7). The level of fractional polarisation 
depends on location within the lobe. The northern outer lobe 
showed higher levels (40-60 %) of FP around the outer parts 
and lower levels (7-20%) towards the centre. The southern 
outer lobe showed similar levels but with a more "patchy" 
structure, perhaps caused by a clumpy composition, consis- 
tent with the 20% depolarisation detected by |Konar et al.| 
( 2006 ) . The inner lobes were slightly less polarised but still 
with the higher levels tending towards the sides of the lobes. 



4 RESULTS - B1834+620 

Similarly, B 1834+620 was detected at each of the six fre- 
quency/resolution combinations although much of the ex- 
tended emission was resolved out in the high-frequency A- 
array images. Images at each frequency are shown in Fig. [3] 
with inset plots overlaying the right-hand panel to show the 
inner lobes at high resolution. All the lobes are extended 
along the axis of the large-scale source and the 1.40-GHz 
images hint at some possible emisson extending from the 
outer lobes back towards the inner lobes. This source has 



previously been studied bylSchoenmakers et al. (2000b) and 



here we extend their analysis using new, higher-resolution 
images. 

We used the VLA images to obtain fluxes, widths 
and lengths of the lobes, again typically using the lower- 
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resolution images, and to create spectral index images 
(Fig. [4]). The values of a appear reasonable and of the or- 
der ~ —1 for lobes and flatter in the core; however, the 
spectrum in the 4.86-GHz/8.46-GHz image appears unex- 
pectedly steep (—3 to —4) in places suggesting that some of 
the 8.46-GHz emission has been resolved out. We therefore 



use the fluxes of Schoenmakers et al. (2000b) in Section 7 



The discrepancy affecting the spectral index image also 
appears to affect our estimates for the rotation measure; we 
obtain a very low value of 0.4 rad m -2 which is inconsistent 
with the higher value of 55-60 obtained by Schoenmakers 
et al.|d2000b] . Given the much greater number of frequen- 
cies used by these authors, and the lower resolution of their 
images, we adopt their value for correction of our images. We 
note that unpublished data from the Westerbork Synthesis 
Radio Telescope further support this higher value of rotation 
measure (de Bruyn et al. in prep.). The resultant fractional 
polarisation image at 1.40-GHz is plotted in Fig. [5] with in- 
set plots showing the inner lobes at 4.86-GHz. Contrary to 
B 1450+333, the electric vectors are approximately perpen- 
dicular to the jet axis in each of the four lobes, with the 
main region of disruption taking place around the leading 
edge of the northern outer lobe. This location corresponds 
with the bright hotspot, an uncommon feature of DDRG 
objects, although common to detect in more conventional 
FRII galaxies. The degree of fractional polarisation is less 
extreme in this object, typically around 20% in all lobes, al- 
though possibly a little higher in parts of the northern outer 
lobe which, again, is likely to be related to the hotspot. 



5 A MODEL FOR THE OUTER LOBES 

The standard model for the evolution of the large-scale 
structure inflated by jets in FRII-type radio galaxies was 
first proposed by Scheuer (1974) and Blandford & Rees 



( 1974 ). The momentum transported along the jet is balanced 
by the thermal pressure and ram pressure of the receding 
medium in the vicinity of the host galaxy. The momentum 
balance leads to the formation of a strong shock at the end 
of the supersonic jet. After passing through the jet shock the 
jet material inflates a lobe or cocoon of low density, but high 
pressure, gas. This material protects the jet flow from tur- 
bulent disruption through contact with the higher density 
ambient medium. 

Falle ( 1991| showed that an initially overpressured jet 



will come into pressure equilibrium with its own lobe. For a 
jet propagating into an ambient medium with a power-law 
density distribution, i.e. 



■(;)" 



(1) 



where r is the distance from the centre of the distribution 
and p and a are constants, the growth of the lobe will be 
self-similar (Kaiser fc Ale xander|1997 ). The size of the lobe 
is then related to the age of the jet flow, t, by a power-law. 

We can also calculate the radio luminosity of the lobe 
of a given size (e.g.|Kaiser et al.|1997" Blundell & Rawlings 



2000 |Manolakou fc Kirk|2002| ). This takes into account the 
cumulative effects of energy losses of the radiating electrons 
due to the adiabatic expansion of the lobe, the emission of 
synchrotron radiation and the inverse Compton scattering 
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Figure 5. Fractional linear polarisation (FP) images of 
B1834+620 using the VLA:B-array data only; FP is represented 
by the greyscale, which ranges from 3cr to the map peak in 
each case. The "wedge" indicates the greyscale in the main im- 
age only. The main image is at 1.40-GHz and the inset plots of 
the inner lobes are at 4.86-GHz. The vectors indicate the direc- 
tion of the electric field, following correction for RM (=55 rad 
m -2 ), and their length indicates the linear polarisation (1 ;/ = 
4.81 x 10 -5 Jy beam -1 ). The outer lobes are strongly polarised, 
typically up to ~ 20%. The electric vectors are perpendicular to 
the source axis in all lobes, with the exception of the hotspot in 
the north outer lobe. 



of Cosmic Microwave Background (CMB) photons. Here we 
use the model of the source dynamics by Kaiser & Alexan- 
der (1997) and that for the radio synchrotron emission of 
the lobe by |Kaiser et al.| ([r997). This combined model was 
recently summarised and discussed in Kaiser & Best (2007) 
(see also Kais er &; Best||2008| ) and we will use the model in 
the form presented there. 

In the radio maps of the two DDRGs there is little or no 
evidence for active hotspots in the outer lobes. This suggests 
that they are no longer supplied with energy by active jets. 
Here we assume that the energy supply to the outer lobes 
only stopped recently in terms of the total source age, t. 
The dynamics of the lobes are then still governed by the 
same equations as during the active jet phase lasting a time 
£j < t. However, no freshly relativist ic accelerated electrons 
enter the lobes after tj + tt, where tt is the time it takes 
the last jet material ejected by the central active galactic 
nucleus (AGN) at tj to reach the end of the still growing 



lobe ( |Kaiser et aL]|2000| |Brocksopp et al.||2007| ). We note 



that the parameter tj did not appear in the standard model, 
referenced above, but is introduced here as a free parameter. 

The combined model involves a number of parameters. 
These can be divided into three groups. Set parameters 
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are fixed to plausible values by considering the underlying 
physics or observations independent of the radio observa- 
tions. The value of observed parameters are derived from 
the observational radio data of individual DDRGs discussed 
in this paper. Free parameters are those derived by making 
the model self-consistent and taking into account all obser- 
vational information. Table [2] summarises all parameters and 
we discuss the more important ones below. 

The adiabatic index of the lobe material depends on 
its composition. For lobes filled only with magnetic fields 
and a relativistic pair plasma, we would expect a relativistic 
equation of state with Ti = 4/3. However, in the following we 
find that the energy density of the lobes of B 1450+333, and 
to a much lesser extent those of B 1834+620, may contain a 
significant contribution from other, non-radiating particles, 
i.e. k ^ (where k is the ratio of the energy densities of non- 
radiating particles and relativistic electrons). Also, if the 
power-law energy distribution of the electrons does indeed 
extend down to Lorentz factors of around unity, then most 
particles will not be highly relativistic and hence we set Ti = 
5/3. The limits of the initial electron energy distribution do 
not strongly influence the results as long as m ~ 2, but 
for steeper power-laws the lower energy cut-off may become 
more important. This is the case for the inner lobes of the 
DDRGs discussed here and we will return to this point in 
Section 16.11 

The way in which the ratio of the energy density of the 
magnetic field in the lobe and that of all particles, R, is de- 



fined in Kaiser et al. ( 1997) and|Kaiser fc Best (2007) implies 



that the relativistic electrons and the magnetic field are not 
fulfilling the minimum energy condition at any point in time 
when k ^ 0. However the magnetised plasma in the lobes of 
FRII-type sources appears to be close to minimum energy 



conditions (e.g. Croston et al. 2005). Hence we enforce min- 
imum energy conditions by setting R — 3/ [4 (k + 1)]. In this 
way, all model equations involving R and k remain valid. 

The outer lobes of the DDRGs discussed in this paper 
are very large. The model therefore predicts the sources to 
be old. Any inclination of the lobes towards our line of sight, 
implying a viewing angle < 90°, would result in projection 
with the actual lobe sizes and their ages increasing even fur- 
ther. However, this is a small effect as long as > 45°. While 
the observational data do not rule out projection effects, we 
assume that the lobes are lying in the plane of the sky. 

The density distribution of the ambient gas is described 
by equation ([!]). The model therefore only depends on the 
combination pa 13 , but not on a and p individually. Here we 
choose to fix a and adopt p as a free parameter. Any change 
suggested by the model for p may be interpreted as a change 
of a instead as long as the combination pa 13 retains the value 
suggested by the model. The exponent of the power-law f3 
is set to 1.5, appropriate for a wide range of gaseous en- 
vironments from elliptical galaxies to galaxy clusters (e.g. 
Fukazawa et al. 2004). We also investigate the possibility 



that the outer radio lobes extend beyond the local environ- 
ment of the host galaxy or cluster into a more uniform Inter 
Galactic Medium (IGM). In this case we set /3 = 0. 

The length of the lobes and the ratio of their length 
and their width, the lobe aspect ratio A, can be determined 
from radio observations with sufficient resolution to resolve 
the lobe perpendicular to the jet axis. However, at high ob- 
serving frequencies, allowing sufficient spatial resolution, the 



lobes are often not detected all the way from the hotspots 
to the core. While we can still determine the length of the 
lobes, it is difficult to measure their width. Most measure- 
ments of A must therefore be viewed as upper limits. 

Measuring the total radio luminosity density at a given 
observing frequency requires low spatial resolution to avoid 
resolving out flux on the largest scales. At the same time we 
require sufficient resolution to distinguish emission from the 
inner and outer lobes. For both sources we find measured 
fluxes at a range of frequencies that avoid both problems. 

The power-law exponent for the initial energy distribu- 
tion of the relativistic electrons may, in principle, be deter- 
mined from the observed spectral slope of the lobes at low 
frequencies where radiative energy losses of the electrons 
have not yet significantly changed the energy distribution. 
This is straightforward in the case of young sources like the 
inner lobes and below we will use the value for m derived 
from their observed spectral slope. A similar determination 
of m is difficult to achieve for older sources like the outer 
lobes of the DDRGs (e.g. |Konar et al.[|2006| ). We show below 
that the observed spectra of the outer lobes are consistent 
with a comparatively flat energy distribution with ra = 2. 

We use the observational data to constrain the values 
of four free parameters, the jet power, Q, the normalisation 
of the density distribution of the ambient medium, p, the 
source age, t, and the duration of the jet flow, tj. Assuming 
that we have n measurements of the radio flux at the same 
number of frequencies spread out over a large range plus 
the length of the lobe, then the degrees of freedom of the 
model are given by n — 3. For the inner lobes we assume the 
jets are still supplying these structures with energy and so 
the number of free parameters is reduced to three, because 
t — tj. Also, we show below that the spectra of the inner 
lobes are consistent with power-laws. Hence the individual 
flux measurements are not independent and they can only 
be counted as a single measurement. For the inner lobes we 
therefore have only two observational constraints, but three 
free parameters, resulting in relations between pairs of free 
parameters instead of specific values for each of them. 



6 MODELS FOR THE INNER LOBES 

The spectra of the inner lobes, up to an observed frequency 
of at least 8.46 GHz, are consistent with a power- law de- 
scription. This implies that the electrons emitting in the 
observed frequency range have not suffered significant radia- 
tive energy losses since they were accelerated to relativistic 
velocities. The inner lobes must therefore be young. This re- 
quirement is consistent with our assumption that the inner 
lobes represent a second activity phase of the jet flow from 
the central AGN after the energy supply to the outer lobes 
shuts down. The inner lobes then cannot be older than a 
few 10 7 years, because for an older age the emission from 
the outer lobes should have faded below the fluxes at the 
higher observed frequencies. Given the observed close align- 
ment of the inner lobes with the outer lobes, we assume that 
the inner lobes are fully contained within the volumes of the 
outer lobes. We propose two different interpretations for the 
existence and observed properties of the inner lobes. 
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Adiabatic index of magnetic field 
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Minimum of initial electron Lorentz factor 


Tmin 


1 


Set 


Maximum of initial electron Lorentz factor 


7max 


10 8 




Viewing angle 


e 


90° 




Exponent of external density distribution 


P 


0, 1.5 




Core radius of external density distribution 


a 


2kpc 




Friprcrv Hprmitv nf thp flMR at 7 — D 


W CMB 


4 x 1 -14 T m -3 

Tfc /\ 1U O 111 




Ratio of energy densities of non-radiating particles and relativistic electrons 


k 


0, 100 




Length of lobe 


D 






Aspect ratio of lobe 


A 




Observed 


Observing frequency 


V 






Radio luminosity density 


L v 






Cosmological redshift 


z 






Exponent of initial power-law energy distribution of relativistic electrons 


m 






Jet power 


Q 




Free 


External density at core radius 


p 






Source age 


t 






Duration of jet flow 


h 





Table 2. Summary of model parameters for the standard FRII model. The notation is that used in Kaiser & Best ( |2007| >. 



6.1 Standard FRII model 

The first interpretation is that the same model used to ex- 
plain the outer lobes also applies to the inner lobes. The jet 
flows inflating the inner lobes must then end in strong shocks 
giving rise to particle acceleration. The observations do not 
show much evidence for luminous radio hotspots which could 
be identified as the jet shocks. However, given our limited 
spatial resolution, we cannot rule out the presence of com- 
pact structures with an enhanced radio luminosity at the 
end of the inner lobes. 

The model is identical to that developed in Section [5] for 
the outer lobes. However, we can make some simplifications. 
We assume that the inner lobes are currently still supplied 
with energy by active jets. Hence tj = t and we reduce the 
number of free parameters by one. Also, the power-law shape 
of the observed spectrum implies insignificant radiative en- 
ergy losses of the relativistic electrons. The inner lobes will 
be well described by the much simplified model in the adi- 
abatic regime as described in |Kaiser" fe Best| (T2007). For an 
observed lobe length of D and radio luminosity density L u 
it is straightforward to derive relations between the source 
age of the inner lobes, ti, and the jet power, Q, 



( 1 + e )4 jD 3(m+l) jL 4 



ftfS 



-. l/(m+5) 



m + 1 (5-/3)(m+l)/3 



(2) 



where e, /l, f P and c\ are constants defined in Kai ser fe| 
Best (2007). We can also find a relation between the scale 



density of the ambient gas, pi, and the jet power, 

^5-/3 

Pi 



(I)' 



a 



-13 



Qtf. 



(3) 



The observed power-law spectra have steep slopes and 
in the absence of significant radiative energy losses these 
steep slopes translate into large values for the exponent of 
the initial energy distribution of the relativistic electrons, m. 
Large m may indicate inefficient particle acceleration and 
they make the model results more dependent on the low en- 
ergy cut-off of the electron energy distribution as most of 



the energy of the electrons is stored in the low energy end of 
the distribution. However, our conclusions below do not de- 
pend on our assumptions for 7 m in- In general, this standard 
FRII interpretation of the inner lobes requires substantially 
higher mass densities within the outer lobes than can be 
explained with the material transported along the jet flows 
that inflated these structures. 



6.2 Bow shock model 

For the alternative model we assume that the density and 
pressure inside the outer lobes is insufficient to significantly 
slow down the advance of the inner, young jets. In this 
case, the inner jets do not develop very strong shocks at 
their ends or innate substantial lobes. We do not expect the 
ends of the jets to be sites of efficient particle acceleration 
and so the inner jets do not directly produce much radio 
emission. However, they push aside and compress the con- 
tents of the outer lobes. The advance of the inner jets is 
fast and therefore leads to the formation of a bow shock 
within the outer lobes. This bow shock can re-energise the 
relativistic electrons and also strengthens the magnetic field 
by compressing it. We identify the radio synchrotron emis- 
sion observed from the inner lobes with that produced by 
the shocked and compressed outer lobe material. The model 



was discussed briefly in Brocksopp et al. (2007), where we 
also point out that bow shock structures, but no lobes di- 
rectly inflated by jets, are found in numerical simulations 



of restarting jets (Clarke & Burns 1991). Here we expand 



the description of the model to see whether our observa- 
tions of B 1450+333 and B 1834+620 are consistent with its 
predictions. A very similar model was recently developed in- 
dependently by Safouris et al. (2008) to describe the inner 



lobes of the DDRG PKS B1545-321. 

In the restframe of the bow shock propagating through 
the outer lobe ahead of the inner jet, the momentum sup- 
plied by the jet must be balanced by the momentum of the 
receding lobe material. We must also include the pressures 
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of the jet material, pj, and of the lobe material, p\, in this 
balance and so (e.g. Bicknell 1994| 

(4) 



(7?b$bWj 4i = (7b/^b^i +Pi) -4b, 

where /3j is the speed of the bow shock in the restframe of 
the outer lobe material and /?jb is the speed of the jet ma- 
terial in the restframe of the bow shock. 7b and 7jb are the 
corresponding Lorentz factors. From the relativist ic trans- 
formation of velocities we get jfbPfb — 7hPf (A — AO 2 - The 
jet has a cross-section of Aj and the bow shock has a surface 
area of Ah. 

The relativistic enthalpy w is defined as w = e + p with 
e the internal energy density, including the restmass energy 
density, and p the pressure. For simplicity we assume that 
all gases can be described by an equation of state of the form 
p = (r — 1) (e — pc 2 ^ , where Y is the ratio of specific heats 
and p the restmass density. It is convenient to introduce the 
ratio 1Z — pc 2 / (w — pc 2 ) . For large 1Z the energy density of 
the gas is dominated by its restmass energy density, while 
for small values of 1Z the gas is relativistic. In Brocksopp 
et al. (2007) we only considered the latter case. We can now 
write 



(Wj + 1) + 1 



Mi 



7b /3b 



(fti + l) + l 



Tj-1 

The energy flux along the jet is 



P\A h . 



(5) 



Q = (t^J - Pjc) facAj. (6) 

Substituting for Wj and p } in terms of pj and 1Z] and solving 
for the jet cross-section gives 



Q 



Ij iy Vi -f ■ R: (-f \) 



(7) 



Note that upon substituting this expression for Aj into equa- 
tion ([5]), the jet pressure terms cancel. 

We now turn to the properties of the material in the 
outer lobe and find an expression for 7Z\. The jet inflated 
this lobe during a time tj. If its volume is V\ and it only 
contains material that was supplied by the jet, i.e. mixing 
across the lobe surface can be neglected, then 



Pi 



(8) 



where M is the restmass transport rate of the jet. M is given 
by 

M = P AcA j . (9) 

Substituting for pj in terms of 1Z } and p } as well as substi- 
tuting for Aj we find 

Q Hi 



M 



Since 



Hi 



I pic 2 



Pi 



(10) 



(11) 



we can now find the value of 1Z \ , if we know p\ and the other 
quantities appearing in equations |5]) and |7| . We obtain p\ , 
Q, and V\ from the model for the outer lobe, while we can 
use the observed geometry of the inner lobes to constrain the 



surface area of the bow shock, Ah- Assuming either a matter 
dominated, 7£j ^> 1, relativistic, 7Zj <C 1, or balanced, 7Zj ~ 
1, jet we can then find the speed of the bow shock, /3b, as a 
function of the Lorentz factor of the jet, 7j , from the implicit 
equation |5]). 

The material in the outer lobe passes through the bow 
shock driven by the inner jet. In this process, the mass, mo- 
mentum and energy of the lobe material must be conserved. 
This leads to the continuity equations across a relativistic 
shock (e.g. |Landau &; Lifshitz 1987), 

Pbjbpl = /3s7s/>s, 



A>7bV 



PbUbWb + Ps 



(12) 



where a subscript 's' denotes quantities describing the 
shocked material behind the bow shock. All velocities are 
measured in the restframe of the bow shock. If we assume 
that the ratio of specific heats of the material in the outer 
lobe, Ti, does not change during the shock passage, then we 
can re-arrange the set of equations above to give 



Pbjb 



and 



7b r t 

7s Ti - 



-A* (A> - A) 7b 7s 



A) 7b - A7g 

1 + 



(13) 



1 1 ;f3bjb (Pb 



V 



■Ps) (1+ ft + 1 



pi- 



(14) 



From the implicit equation (13) we obtain the velocity of 
the shocked gas with respect to the bow shock, /3 B . Substi- 



tuting this into equation (14) we arrive at the pressure of 



the shocked gas p s which may be compared to the pressure 
estimated from the radio synchrotron emission of the inner 
lobes by using the minimum energy condition. 



7 APPLICATION OF THE MODELS 

We now apply the models detailed above to the observa- 
tional data on B 1450+333 and B 1834+620. 



7.1 B 1450+333 

The radio observations used to constrain the models for the 
inner and outer lobes of B 1450+333 are summarised in Ta- 
ble [3] We follow |Konar et ah] §006 ) and assume a typical 
error for the flux measurements of 7% at all frequencies, 
except at 0.240 GHz where the error is assumed to be 15%. 



7.1.1 Outer lobes 

For the outer lobes we have flux measurements at six fre- 
quencies for each lobe. According to the discussion in Sec- 
tion [5] this implies that we have three degrees of freedom 
in the model fitting. We investigate three individual mod- 
els. The first assumes a power-law density distribution the 
jets expand into with /3 = 1.5 and no energy in any non- 
radiating particles in the lobes, i.e. k = 0. We then allow 
for non-radiating particles in the lobe by setting k = 100. 
Finally we assume that the outer radio lobes extend beyond 
the confines of the atmosphere of the host galaxy into a uni- 
form intergalactic medium with f3 = 0. For this model we 
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53 


43 


55 


48 


52 


43 


tj /Myr 


48 


34 


50 


39 


46 


33 


p/10- 24 kgm- 3 


3.0 


0.58 


1.3 x 10- 3 


2.2 x 10- 4 


170 


34 


x 2 


0.68 


0.86 


0.35 


0.73 


0.70 


0.75 


pi/10- 15 Jm- 3 


3.1 


1.5 


3.5 


1.7 


180 


88 


Si/nT 


0.071 


0.049 


0.075 
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0.071 


0.049 
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6.4 x 10- 7 


8.2 x 10- 7 


5.4 x 10- 7 


7.2 x 10- 7 


3.7 x 10- 5 


4.8 x 10- 5 


Pi/10- 24 kgm- 3 


0.023 


2.0 x 10- 4 


0.036 


2.0 x 10- 4 


7.4 x 10- 4 


3.6 x lO -6 


ti/Myr 


2.4 


0.41 


2.9 


0.54 


0.19 
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Table 4. Model parameters for the outer lobes of B 1450+333 for the model spectra with the best fit. The model parameters for the 
inner lobes for the standard FRII model, t{ and pi, evolving inside the respective outer lobes are also given. 



B1450+333: 


North 


South 




Outer 


Inner 


Outer 


Inner 


Lobe length, D / kpc 


635 a 


79.8 b 


781 a 


83. 7 b 


Aspect ratio, A 


1.8° 


3.9 d 


3 C 


5 d 


Flux densities / mJy 


0.240 GHz 


1055 


123 


602 


48 


0.334 GHz 


902 


96 


563 


47 


0.605 GHz 


606 


69 


370 


29 


1.287 GHz 


272 


37 


164 


17 


1.365 GHz 




38 




14 


1.40 GHz 


258 




145 




4.86 GHz 


71 


16 


32 


4.9 


8.46 GHz 




9.5 




2.6 



Table 3. Observed properties of the four lobes in the DDRG 
B 1450+333 used in the modelling. Additional flux densities have 

|2006| . 



been taken from the compilation by 



Konar et al. 



a VLA in B-array at 1.40 GHz 

b VLA in A-array at 4. 86 GHz 

c GMRT at 0.605 GHz ( |Konar et ah] [2006} 
d VLA in A-array at 1.40 GHz 



implicitly assume that the lobes have extended beyond the 
atmosphere of their host galaxy for a significant fraction of 
their age. Table El shows the fitting results for these three 
models and Fig. [6] shows the model spectrum for the case 
f3 = 1.5 and k = 100 compared to the observations. 

In all three cases we get very good agreement of the 
model spectra with the observational data, despite the fact 
that the observed spectral slope at high frequencies is con- 
siderably steeper than the a = 0.5 expected from a model 
with m = 2 in the absence of radiative energy losses of the 
electrons. The fact that x 2 < 1 for all three cases suggests 
that our assumed errors of the flux measurements, 15% at 
0.240 GHz and 7% at the other frequencies, are too conserva- 
tive. Also, based on the % 2 values alone, none of the models 
is formally preferred over the others. 

While we cannot formally reject any of the models, a 
closer inspection of the model results shows that both mod- 
els with k = require a very low gas density in the source 
environment. We can convert the gas density measured from 




l.o 

Frequency / GHz 



10.0 



Figure 6. Comparison of model spectra with observations for 
B 1450+333 in the source restframe. Circles show the data for 
the northern lobes, squares are for the southern lobes. The upper 
data points and lines are associated with the outer lobes while 
the lower data points and lines show the inner lobes. The solid 
lines show the model predictions for the outer lobes, the dashed 
lines show the power-law fits to the spectra of the inner lobes. 



to our choice of /3 = 1.5 and a = 2 kpc. The lowest density in 
this sample is that of NGC 3607 with p-8x 10~ 23 kgm -3 
followed by NGC 720 with p ~ 1.40 x 10" 22 kg m - 3 . The 
densities required by our model with /3 = 1.5 are at least a 
factor 25 lower. Also, for /3 = 0, when we assume that the 
outer lobes extend beyond the gaseous atmosphere of the 
host galaxy, the inferred density is at least a factor 14 be- 
low the critical density of the universe at this redshift. We 
consider it therefore unlikely that the models with k = 
accurately describe the outer lobes of B 1450+333. 

We experimented with changing the values of the vari- 
ous set model parameters to arrive at a more realistic density 
for the source environment. The only model parameter forc- 
ing a significant change of p is k, the ratio of the energy 
stored in non-radiating particles inside the outer lobes. Set- 
ting k = 100, we found an acceptable fit with a density for 
the host environment consistent with the findings of | Jetha] 
et al. (2007) for galaxy groups. The age of the outer lobes 



X-ray observations of galaxy groups by Jetha et al. (2007) 



does not change much for this model because the observed 
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spectral shape requires significant ageing of the radiating 
electrons. However, the larger energy requirements force a 
large increase in the jet power, Q. While still larger values 
for k are conceivable, given that for k = 100 the density 
inferred by the model is still at the lower end of the range 



100 



of the sample of |Jetha et al. (2007), the required increase in 
Q would put the jet power of B 1450+333 beyond the maxi- 
mum inferred for radio galaxies of type FRII (e.g. Rawlings 
|&; Saunders[ 1991). Allowing m to be a free parameter can 
help to reduce the high required values of Q (compare with 
e.g. Machalski et al. 2007), although at risk of making the 



model too flexible. However, increasing m too much results 
in the outer lobes becoming so young that there is insuffi- 
cient time for the inner lobes to develop. 

The strength of the magnetic field inferred by the model 
is independent of the value of k while the pressure in the 
lobes increases substantially for larger k. This is just a reflec- 
tion of the assumption, implicit in our model, of conditions 
close to equipartition inside the lobes. A given radio lumi- 
nosity originating in a given lobe volume fully determines 



the strength of the magnetic field (e.g. Longair 1994). The 
magnetic field and the synchrotron emitting electrons both 
contribute to the lobe pressure, p\, but for k = 100, p\ is 
dominated by the energy density of the non-radiating par- 
ticles. 

The uncertainties in the flux measurements at the fre- 
quencies used to fit the model spectra give rise to uncer- 
tainties in the free model parameters. The density in the 
host galaxy environment is associated with the largest un- 
certainty. Fig. [7] shows the confidence contours in the jet 
power - density model plane. The contours for the outer 
lobes overlap considerably with respect to the range in jet 
power preferred by the model, but the model indicates a 
higher density to the north of the host galaxy. Note that in 
this and the following figures the confidence contours would 
cover smaller ranges for the smaller errors on the flux mea- 
surements suggested above. The quality of the data makes 
it particularly difficult to constrain the source parameters; 
however, for further discussion of the discrepancies between 
model parameters for each pair of lobes, see Machalski et al.| 
( 2009] ). 

Fig. [8] shows the confidence contours in the jet power - 
source age plane. Both Q and t should be the same for both 
lobes and it is reassuring that the overlap of the confidence 
contours in this plane is large. The model constrains the 
source age and the jet power roughly to within a factor three. 

Finally, the confidence contours in the jet power - mag- 
netic field strength plane are shown in Fig. [9] As men- 
tioned above, the model strongly constrains the strength 
of the magnetic field in the lobes. This demonstrates that 
the strength of the magnetic field and hence the pressures in 
the lobes are the best constrained parameters of the model. 
The values of B\ and p\ do not change much even for the 
assumption of ft = or for changes of the other set model 
parameters. 



7.1.2 Inner lobes 

The observed radio spectra of the inner lobes of B 1450+333 
are well described by the power-laws 




Jet power / 10' 



Figure 7. Confidence contours for the free model parameters jet 
power, Q, and environment density, p, for the model of the outer 
lobes of B 1450+333 with (3 = 1.5 and k = 100. Upper contours 
for the north outer lobe, lower contours for the south outer lobe. 
The contours from the outside inwards are for confidence levels 
99%, 95%, 90%, 75% and 50%. Black dots show the location of 
the best fit parameters summarised in Table ^] 
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Figure 8. Same as Fig. [7] but for jet power, Q, and source age, 
t. The contours more extended towards lower jet powers show 
the model parameters for the north outer lobe. The other set of 
contours shows the south outer lobe. 
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for the northern inner lobe and 

( v Y 



: 5.8 I 



Jy 



(15) 



(16) 



VMHz/ 

for the southern inner lobe. Fig.|6]shows these fits translated 
to the source restframe. We can now apply the modified 
standard FRII model described in section l6Tl to the obser- 
vations by setting m — 2.42 for the northern inner lobe and 
m — 2.68 for the southern inner lobe. This assumes that 
the radiative energy losses of the relativistic electrons in the 
inner lobes are not significant enough to affect the spectra. 

As discussed in section |6J") the set of model parameters 
Q, pi and U cannot be uniquely fitted to the observations 
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Figure 9. Same as Fig. [TJ but for jet power, Q, and magnetic field 
strength in the lobe, B\. The upper contours show the parameters 
for the north outer lobe, the lower contours for the south outer 
lobe. 



in this case. Hence we take the value of Q from each of the 
models fitted to the data of the respective outer lobe and 
calculate pi and U for this value of Q from equations |2]) 
and |3j). Here, we make the implicit assumption that the 
power of the jets that innate the inner lobes is equal to the 
power of the jets that produced the outer lobes, although it 
is indeed possible that this assumption is not valid (see e.g. 



Schoenmakers et al. 2000a Machalski et al. , 2010). We also 



assume that the density distribution inside the outer lobes, 
which provide the environment the inner lobes are expand- 
ing into, is uniform, i.e. f3 = 0. The results are summarised 
in Table 

Clearly, the standard FRII model requires young ages 
for the inner lobes, t\. All values of t\ are consistent with the 
requirement that the inner lobes must have been inflated 
since the jet flows to the outer lobes stopped, i.e. in every 
case t[ < t — ty The young ages also require the advance 
speed of the leading edge of the inner lobes to be fast. For the 
models of the outer lobes with k — 0, these advance speeds 
averaged over the age of the inner lobes lie between 0.1c 
and 0.7 c. However, for the preferred models with k = 100 
the average advance speeds become superluminal. Such an 
unphysical situation can only be avoided if the power of the 
jets inflating the inner lobes drops compared to the jet power 
associated with the outer lobes. 

Another problem for the standard FRII model arises 
from the densities required by the model for the inner lobes. 
Assuming a jet power Q, the model makes a prediction for 
the density of the gas surrounding the inner lobes, p x . Since 
the inner lobes expand inside the outer lobes, this density 
Pi must be equal to the density inside the outer lobes. Com- 
bining equations ([8} and ( 10 ) we can derive an estimate for 



the density inside the outer lobes, p\. The density p\ arises 
solely from the material that was transported along the jets 
that inflated the outer lobes. It does not take into account a 
possible 'contamination' of the outer lobes by mixing in ma- 
terial across the outer lobe boundaries. For simplicity Table 
[4] lists the values for p\ (£ = 1) after setting the factor 



(17) 



equal to unity. For a self-consistent model without mixing 
across the outer lobe boundary we require pi — p\. It is 
straightforward to show that this results in 



7j 



i-7i 2 + M£ = i)/pi' 



(18) 



By definition Tt } > and so we find a limit on the bulk 
Lorentz factor of the jet 



7j 



< 



Pi (£ = i) 



+ 1. 



(19) 



The models with k — predict p\ (£ = 1) <C Pi- While 
for models with k = 100 the two density estimates are more 
comparable, we argued above that in these cases we expect 
a reduced power for the jets inflating the inner lobes which 
implies a higher value for pi more comparable with that for 
the models with k — 0. Therefore all models predict rather 
slow jets with y } ~ 1, dominated by the transport of kinetic 
energy. Such slow jet velocities are unlikely and so to make 
the standard FRII model work we need to invoke significant 
mixing of high density gas from the source environment into 
the outer lobes. This is of course consistent with the find- 



ings of |Kaiser et al. ( 2000| ), who applied the same model to 
DDRGs. 

We now turn to the bow shock model for the inner lobes. 
To calculate the pressure inside the inner lobes, p s , as pre- 
dicted by the dynamical model described in Section [6. 2 1 we 
use the model for the outer lobes with f3 — 1.5 and k — 100. 
This sets the jet power, Q, the volume of the outer lobes, V\, 
the duration of the jet flow inflating the outer lobes, tj, and 
the pressure inside the outer lobes, pi. From the observed 
geometry of the inner lobes we determine the area of the 
leading surface of the bow shock, Ah- We assume a simple 
cylindrical geometry so that Ah = k[D/ (2A)] 2 . For given 
values of 1Z } we can now calculate predictions for the pres- 
sure inside the inner lobes, p s , as a function of the Lorentz 
factor of the jet, j y 

Fig. [l0| shows the results for three values of 7Z y Rela- 
tivists jets dominated by the internal energy of the jet mate- 
rial, i.e. 7Zj = 0, predict the highest pressure inside the inner 
lobes for a given jet Lorentz factor, while the cold, matter 
dominated jets with IZj — >• 00 result in the lowest pressures. 
The endpoints of the curves at small 7j are determined by 
the requirement that the shocked material behind the bow 
shock cannot move faster in the lobe restframe than the jet 
material itself. At large 7j the curves terminate where equa- 
tion ( 14 ) implies that there is no bow shock formed ahead 
of the jets, i.e. p s = p\. 

For a fixed jet power, Q, and value of 7Zj the jet cross 
section, Aj, is a decreasing function of 7j (see equation [7]) . 
In other words, a faster jet achieves the same energy trans- 
port rate through a smaller cross section. Hence, equation 
|5]) predicts a slower advance velocity for the bow shock, /3b, 
for increasing 7j. This in turn reduces the compression ra- 
tio of the bow shock and the decreasing predicted pressure 
inside the inner lobes seen in Fig. |10| Also, a slower bow 
shock implies a slower velocity for the shocked material be- 
hind the bow shock ahead of the jet. At the endpoints of 
the lines for large 7j, the shocked material is almost at rest 
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B1834+620: 



North 



South 




Figure 10. The pressure in the inner lobes of B 1450+333 pre- 
dicted by the bow shock model as a function of the Lorentz factor 
of the jet. Light lines show the results for the north inner lobe, 
heavy lines for the south inner lobe. The dashed lines show re- 
sults for relativistic jets dominated by the internal energy of the 
jet material, i.e. 7Zj = 0, while dash-dotted lines show rest mass 
dominated jets with 7Zj — > oo. Solid lines show the results for bal- 
anced jets, 7£j = 1. The horizontal dotted lines show the pressure 
inside the inner lobes derived from observed radio flux assuming a 
cylindrical geometry of the emission region and minimum energy 
conditions and setting k = 10. 



inside the lobe restframe. Accordingly, the bulk velocity of 
the jet relative to the shocked lobe material increases from 
left to right in Fig. |10| It is unlikely that the jet material 
can decelerate without forming a strong shock when this 
relative velocity is large. The absence of observed hotspots 
in the inner lobes and the problems of the standard FRII 
model discussed above therefore argue for solutions towards 
lower values of 7j . 

Assuming minimum energy conditions inside the in- 
ner lobes, we can derive an independent estimate for the 
strength of the magnetic field (e.g. Longair 1994). We as- 
sume a cylindrical geometry of the inner lobes to calculate 



their volume. We have argued above that the radio spectrum 



is well fitted by the power-laws given in equations (15) and 



(16). The estimate of the minimum energy magnetic field, 
£>min, should therefore not be significantly affected by en- 
ergy losses of the radiating electrons and we use the spectral 
slope of the fitted power-laws. We derive two values for Bmin 
for each of the inner lobes. The 'conservative' value assumes 
that there is no radio emission outside the redshift-corrected 
frequency band across which we have observations of the in- 
ner lobes (see Table |3j). The 'traditional' value assumes that 
the radio spectrum extends between the essentially arbitrary 
limits 10 MHz to 100 GHz. The traditional value will always 
exceed the conservative one. In practice the conservative and 
traditional values for .Bmin are not very different. For the in- 
ner north lobe we find B m i n equal to 0.22 nT or 0.32 nT, 
respectively. For the inner south lobe we get 0.18 nT and 
0.28 nT. 

Fig. [lO] also shows the pressures inside the inner lobes 
derived from the traditional values of B m i n under the as- 
sumption that the non-radiating particles in the inner lobes 
store ten times more energy than the relativistic electrons, 





Outer 


Inner 


Outer 


Inner 


Lobe length, D / kpc 


713 a 


202 b 


736 a 


194 b 


Aspect ratio, A 


12 a 


14 b 


13 a 


16 b 


Flux densities / mJy 


0.612 GHz 


573 


309 


357 


422 


0.845 GHz 


435 


228 


259 


298 


1.40 GHz 


275 


144 


177 


198 


4.85 GHz 




30.9 




51 


8.46 GHz 


53 


19.4 


32.3 


29.1 



Table 5. Observed properties of the four lobes in the DDRG 
B 1834+620 used in the modelling. All flux densities are taken 
from the compilation by Schoenmakers et al. (2000b). We do not 
use the flux density measurements at 1.395 GHz and 1.435 GHz 
quoted in that paper as these observations are suspected to re- 
solve out some of the flux. 



a VLA in B-array at 8.5 GHz 
b VLA in A-array at 4.8 GHz 



i.e. k = 10. This is lower than in the outer lobes with 
k = 100 which the bow shocks are expanding into. How- 
ever, the pressure derived assuming minimum energy con- 
ditions in the inner lobes exceeds any pressure achievable 
by the bow shock. Reducing the value of k across the bow 
shock implies that the shock channels energy preferentially 
to the relativistic electrons and the magnetic field. Some re- 
acceleration of electrons at the bow shock seems plausible, 
particularly since the observed spectral slopes of the inner 
lobes are different from those of the outer lobes, and so a 
small reduction of k may occur. Note that k must decrease 
more for higher 7j. Unless k is reduced significantly, this 
again argues for moderate values of 7j . 

The bow shock model can explain the observations of 
the inner lobes without the need to invoke significant mixing 
of the outer lobes with the surrounding heavier gas. To avoid 
the formation of a strong termination shock at the end of 
the jets, the model predicts moderate bulk jet velocities with 
7j ^2—^3. The bow shocks driven by these jets then prop- 
agate through the outer lobes with velocities comparable to 
the bulk jet velocities. Despite their considerable size, the 
inner lobes are young in this model with ages of only a few 
10 5 years consistent with the absence of any spectral ageing 
signature from the observed spectra. These young ages are 
also in agreement with the time available for the formation 
of the inner lobes, t — tj , as predicted by the model of the 
outer lobes. 



7.2 B 1834+620 

Tabled summarises the observational data used to constrain 
the models for the inner and outer lobes of B 1834+620. As 
for B 1450+333, we assume a typical error of 7% on all flux 
measurements. 



7.2.1 Outer lobes 

For the outer lobes we have flux measurements at four ob- 
serving frequencies. It follows from Section [5] that there is 
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10 24 : 



Frequency / GHz 

Figure 11. Comparison of model spectra with observations for 
B 1834+620 in the source restframe. Circles show the data for 
the northern lobes, squares are for the southern lobes. The upper 
data points and lines are associated with the outer lobes while 
the lower data points and lines show the inner lobes. The solid 
lines show the model predictions for the outer lobes, the dashed 
lines show the power-law fits to the spectra of the inner lobes. 
The luminosities of the inner lobes have been divided by a factor 
3 for clarity. 




7 10 

Jet power / 10 38 W 

Figure 12. Confidence contours for the free model parameters jet 
power, Q, and environment density, p, for the model of the outer 
lobes of B 1834+620 with (3 = 1.5 and k = 0. Upper contours for 
the north outer lobe, lower contours for the south outer lobe. The 
contours from the outside inwards are for confidence levels 99%, 
95%, 90% and 75%. Black dots show the location of the best fit 
parameters summarised in Table [6] 



only one degree of freedom in the model fit. Again we at- 
tempt fits with three models differing in their assumptions 
about the slope of the external density profile, /3, and the 
ratio of the energies stored in non-radiating particles and 
relativistic electrons, k. We tried models with f3 = 1.5 and 
k = 0, P = and k = and p = 1.5 and k 10. The fitting 
results are summarised in Table [6] and the model spectrum 
for p — 1.5 and k = is compared to the observations in 
Fig. [TTJ 

Both models with f3 — 1.5 achieve good agreement with 
the observational data. The model with f3 = provides a 
significantly weaker fit. The % 2 values for the models with 
P = 1.5 are again below unity and may indicate that the 
error of 7% assumed for the flux measurements may be too 
large. 

All fits use a power-law exponent m — 2 for the en- 
ergy distribution of the electrons when they are injected 
into the lobes. The observed spectral slope is steeper than 
a — 0.5 and so significant ageing affects the spectrum. How- 
ever, the observed spectrum shows very little curvature and 
the flux measurements for the outer lobes can be adequately 
fitted with simple power-laws. In fact, the spectra can be fit- 
ted accurately by the standard FRII model for parameters 
avoiding much spectral ageing, if we assume that m > 2. 
For example, we obtain an excellent fit, \ 2 — 0.38, for the 
spectrum of the north outer lobe by setting m = 2.5. The 
problem with this solution is the very young age predicted 
by the model of t = 6Myr. The outer lobes would have to 
expand with an average velocity of almost 0.4 c and this re- 
quires an extremely high jet power in excess of 10 40 W. It is 
therefore unlikely that such a solution accurately describes 
the outer lobes of B 1834+620. 

The model with p = 1.5 and k = requires a density 
distribution in the source environment consistent with the 
range of properties in the Jetha et al. ( 2007| ) galaxy group 




7 10 

Jet power / 10 38 W 

Figure 13. Same as Fig.^] but for jet power, Q, and source age, 
t. The upper contours show the model parameters for the north 
outer lobe. The lower set of contours shows the south outer lobe. 



sample. In contrast to the situation for B 1450+333, there 
is no need to increase k in order to allow for denser, more 
realistic source environments. The source density required 
for the model with k = 10 is about two times higher than 
the upper end of the range of the sample of | Jetha et al.| 
(2007). Much higher values for k, comparable to the k — 100 
for B 1450+333, would again result in unacceptably high 
jet powers and densities for the source environment. These 
results strongly suggest that k is not very large in the outer 
lobes of B 1834+620. 

Figs. |12| |13| and ^] show the confidence contours on 
the model parameters for the assumed error on the flux 
measurements. Again, the smaller errors argued for above 
would decrease the range of the confidence contours further. 
These plots are analogous to Figs. [7| [8] and [9] for B 1450+333. 
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fc = 0, 


(3 = 1.5 


fc = 0, 


(3 = 


k = 10, 


(3 = 1.5 




North 


South 


North 


South 


North 


South 


Q/10 37 W 


89 


78 


72 


63 


620 


540 


t/Myr 


28 


25 


28 


25 


26 


24 


tj /Myr 


26 


23 


26 


23 


24 


21 


p/10- 24 kgm- 3 


210 


130 


0.075 


0.045 


1200 


770 


x 2 


0.62 


0.95 


2.3 


3.0 


0.49 


0.90 


pi/10- 15 Jm- 3 


18 


13 


21 


15 


120 


84 


Si/nT 


0.17 


0.14 


0.18 


0.16 


0.17 


0.14 


py (e = l)/10- 24 kgm- 3 


1.4 x 10- 4 


1.2 x 10- 4 


1.1 x 10- 4 


9.3 x 10- 5 


9.0 x 10- 4 


7.3 x 10- 4 


Pi/10- 24 kgm- 3 


0.066 


0.15 


0.10 


0.22 


0.012 


0.026 


ti/Myr 


2.8 


3.3 


3.5 


4.1 


0.84 


0.97 



Table 6. Model parameters for the outer lobes of B 1834+620 for the model spectra with the best fit. The model parameters for the 
inner lobes for the standard FRII model, t{ and pi, evolving inside the respective outer lobes are also given. For this source the values of 
ti and pi are calculated for a truncated radio spectrum. See text for a detailed discussion. 
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0.13- 
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Jet power / 10 38 W 

Figure 14. Same as Fig. [12] but for jet power, Q, and magnetic 
field strength in the lobe, B\. The upper contours show the pa- 
rameters for the north outer lobe, the lower contours for the south 
outer lobe. 



The model parameters are constrained to a similar degree 
for B 1834+620. However, in Fig.[l3]the contours for the two 
outer lobes do not overlap. We would expect that the plotted 
parameters, the jet power and the source age, should be the 
same for both lobes on either side of the source. While this 
discrepancy is certainly visible in the figure, it is not very 
large. It is impossible to trace the outer lobes of B 1834+620 
all the way back to the core in our radio maps. Hence our 
possible overestimate of the aspect ratio of one of the lobes 
discussed in section [5] may cause this problem. 

An important constraint arising from the model of the 
outer lobes is the very short time available for the inner 
lobes to expand. They can only be inflated after the jet flow 
to the outer lobes ceases at tj. Hence the age of the inner 
lobes must be less than t — tj ~ 2 Myr. Given the large size 
of the inner lobes of B 1834+620, this implies a very fast 
expansion. We will see in the next section that this tight 
constraint makes the case for the bow shock model for the 
inner lobes even stronger than in the case of B 1450+333. 



7.2.2 Inner lobes 

The observed spectra of the inner lobes of B 1834+620 are 
steeper than any of the other spectra discussed in this paper. 
The spectrum of the north inner lobe is well fitted by 

^ = 350 (ra)" 109jy ' ™ 

while the spectrum of the south inner lobe is close to 

^ = 300 (ra)" 102jy - ^ 

Because of its natter spectral slope, the southern lobe is 
brighter in the observed range. This can be seen in Fig. [TT] 
where the fits are compared to the data in the source rest- 
frame. Note that the luminosities of the inner lobes are 
scaled by a factor 3 in this figure for clarity. 

We apply the standard FRII model to the inner lobes 
of B 1834+620 in the same way as in the case of B 1450+333 
with one important difference. From the power-laws fitted 
to the spectra we find that we need to set m — 3.18 for 
the north inner lobe and m — 3.04 for the south inner lobe. 
Hence the bulk of the energy in electrons in the model is 
stored in those electrons with the lowest Lorentz factors 
around 7min- So far we have used 7min = 1 and so those 
electrons storing the bulk of the energy do not contribute 
to the radio emission. As a consequence the standard FRII 
model predicts ages in excess of 10 8 years for the inner lobes, 
because the jets need to deliver a vast amount of energy to 
the inner lobes. The problem is solved by increasing the low 
energy cut-off of the electron energy distribution such that 
7min = 1000. The introduction of such a high cut-off leads 
to the truncation of the model radio spectrum. In fact, we 
cannot set a value for 7 m in much larger than the 1000 used 
here or this cut-off will start affecting the model spectrum 
in the observed range. We did not introduce an equivalent 
high cut-off in the models for B 1450+333 as for the natter 
spectra of this source such a change makes little difference 
to the modelling results. Also, a further decrease in the age 
of the inner lobes of B 1450+333 exacerbates the problems 
noted above. 

The results of the modelling are given in Table [6] For 
this object the ages predicted for the inner lobes are too old 
for k = to be consistent with the properties of the outer 
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Figure 15. Same as Fig. [To] but for B 1834+620. Note that the 
identical horizontal dotted lines here are calculated for k = 0. 



lobes, despite the introduction of a large value for 7 m in- In 
all cases t- x > t — £j and so for k = we cannot build a 
self-consistent model combining standard FRII models for 
both the outer and the inner lobes. For k = 10 the ages 
predicted for the inner lobes are smaller than t — t- y However, 
in this case the average expansion speed of the inner lobes 
approaches significant fractions of the speed of light. It is 
therefore unlikely that the jets inflating the inner lobes will 
terminate in strong shocks as required by the standard FRII 
model. 

As in the case of B 1450+333 the densities required in- 
side the outer lobes to explain the inner lobes are at least 
an order of magnitude higher than the model for the outer 
lobes predicts. Unless there is significant contamination of 
the outer lobe with dense gas from the source environment, 



equation (19) then implies that the bulk speed of the jet 
must be very low. We conclude that the standard FRII 
model for the inner lobes of B 1834+620 can only work, 
if extensive mixing of gas across the lobe surface has oc- 
curred. This is consistent with our findings for B 1450+333 



and Kaiser et al. (2000) 



Next, we apply the bow shock model to the inner lobes 
of B 1834+620 in the same way as in the case of B 1450+333, 
but for the model of the outer lobes with /3 = 1.5 and k = 
0. The pressure inside the inner lobes as predicted by the 
bow shock model are shown in Fig. |15[ There are several 
differences to the equivalent plot for B 1450+333 shown in 
Fig. [lO] The curves for individual lobes, but for different 
values of IZj are much closer to each other. The curves for 
the north inner lobe are also closer to the south inner lobe. 
This means that for a given jet speed the predicted pressures 
inside the lobes cover a small range independent of the form 
in which the jets transport the energy. For most values of 7j 
the balanced jets with = 1 lead to the highest pressure. 

We also calculate the strength of the magnetic field cor- 
responding to minimum energy conditions in the inner lobes. 
We find a conservative value of 0.53 nT and a traditional 
value of 1.2 nT for B m i n in the north inner lobe. The cor- 
responding values for the south inner lobe are 0.58 nT and 
1.2 nT. The pressure in the inner lobes associated with the 
traditional values of B m i n are shown as the dotted horizon- 



tal line in Fig. |15| Here we have assumed that k = and so 
the horizontal lines represent a firm lower limit for p s . This 
is therefore a much stronger constraint than in the case of 
B 1450+333. In this model we can rule out relativistic jets 
with IZj = 0, because they do not allow high enough pres- 
sures inside the inner lobes. Also, the jet bulk Lorentz factor 
cannot be much greater than 2. The implied age of the inner 
lobes of roughly 0.8 Myr is consistent with the available time 
t—tj rsj 2 Myr. Finally, the inner lobes are not consistent with 
large values of k. As with the outer lobes of B 1834+620, we 
find that the inner lobes cannot contain an energetically very 
significant population of non-radiating particles. 



8 SUMMARY AND OUTLOOK 

We find that the observed properties of the outer lobes of 
both, B 1450+333 and B 1834+620, can be modelled satis- 
factorily with the standard FRII model. The absence of hot 
spots from the radio maps is consistent with our interpre- 
tation that the jets in these sources have ceased to supply 
the outer lobes with energy. The steep radio spectra of the 
outer lobes are consistent with model solutions where the 
acceleration of relativistic electrons stopped a few Myr for 
B 1834+333 to several Myr for B 1450+333 ago. 

Unless the outer lobes of B 1450+333 are embedded in 
an environment with an exceptionally low gas density, we 
find that they must contain a population of non-radiating 
particles that store a significant amount of the energy orig- 
inally transported by the jets. For a ratio of energy stored 
in non-radiating particles to relativistic electrons, k, equal 
to 100 we find a model fit with a realistic gas density in 
the environment. In contrast for B 1834+620 such a large 
value of k would lead to an unrealistically high jet power 
and short source age. While setting k = 10 for this source 
results in an acceptable model fit, the observed properties 
are also consistent with k = 0, a negligible contribution to 
the total energy in the lobe by non-radiating particles. 

In both sources the time available for the inner lobes to 
be inflated, denned as the difference between the total source 
age, t, and the duration of the first jet activity creating the 
outer lobes, tj, is comparatively short. This, together with 
their considerable size, implies that the inner lobes must 
propagate rapidly through the outer lobes, independent of 
the specific model we apply to the inner lobes. 

In both sources the standard FRII model struggles to 
accommodate the observed properties of the inner lobes. The 
model is based on the formation of a strong shock at the 
end of the jets caused by ram pressure balance with the sur- 
rounding medium. The gas density inside the outer lobes 
is simply not high enough to effectively stop the jets un- 
less their bulk velocity is very slow. The outer lobes may 
be 'contaminated' by denser gas from their surroundings, 
but the processes leading to such a contamination are slow 
(Kaiser et al.l 2000). An alternative solution to make the 



standard FRII model consistent with the observations is to 
invoke a significant drop of the power of the jets inflating 
the inner lobes compared to the jets responsible for the outer 
lobes of the same source. While we cannot rule this out for 
B 1450+333, it cannot explain the inner lobes of B 1834+620 
as their predicted age would increase and become even less 
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consistent with the constraint provided by the outer lobes 
of this source. 

The bow shock model for the inner lobes interprets 
them as arising from the emission of relativistic electrons 
compressed and re-accelerated by the bow shock in front of 
the jets inside the outer lobes. In this model the jets in the 
inner lobes do not decelerate significantly, they do not form 
radio hot spots - the presence of hotspots in a DDRG would 
rule out the bowshock model, providing an observational test 
- and the lobes are expected to expand fast. The model is 
consistent with the inner lobes in both sources. Although we 
cannot determine the bulk jet velocity from the models, our 
results suggest that the jet bulk Lorentz factors should be 
around 2. The constraints are much stronger for B 1834+620 
because of its larger inner lobes. In this source we can also 
rule out jets whose energy transport is dominated by the 
internal energy of a highly relativistic jet material. 

In the case of B 1450+333 the bow shock model suggests 
a significant contribution of non-radiating particles to the 
total energy stored in the inner lobes. However, the value 
of k is likely somewhat lower than for the outer lobes in 
this object. For B 1834+620 the inner lobes as well as the 
outer lobes cannot contain non-radiating particles making a 
significant contribution to the total energy. 

A combination of the standard FRII model for the outer 
lobes with the bow shock model for the inner lobes provides 
a self-consistent interpretation for both sources studied in 
this work. With this combined model we do not need to 
invoke a drop in the jet powers or a significant contamination 
of the outer lobes with dense gas as proposed in K aiser et al.| 
(2000). While we cannot rule out these processes, we can 



terruption occurs commonly in radio galaxies, then it should 
be possible to identify DDRG-like structures also in smaller 



explain the observations without them. 

The bow shock model for the inner lobes predicts a very 
rapid propagation of the bow shock through the outer lobes. 
In both sources the working surface of the jets driving the 
bow shocks ahead of them should reach the tip of the outer 
lobes only about 2 to 3 Myrs after they started their ex- 
pansion at the source centres. At this point the source will 
develop new hotspots at the ends of the outer lobes and have 
the appearance of a 'normal', although large, FRII- type ra- 
dio galaxy. This may imply that many, if not all, large radio 
galaxies go through short-lived DDRG phases. We may even 
speculate about a general process acting in radio galaxies 
whereby the jet production is infrequently interrupted in the 
central AGN for short periods of time, say of order 1 Myr. 
The jet channels in the existing radio lobes collapse quickly 
(~ 10 4 years, comparable with the sound-crossing time for 
the jet radius) and after the jets restart, they drive bow 
shocks through the old lobes. During this phase the source 
would be classified as a DDRG. After the bow shocks reach 
the tips of the old lobes, the source becomes a standard 
FRII-type radio galaxy once more. Due to the fast prop- 
agation of the bow shocks the DDRG phase is short-lived 
compared to the total source age. Also, it lasts longer in 
larger sources simply because it takes the bow shocks longer 
to propagate to the ends of already large lobes. This may 
explain why DDRGs are rare and are typically associated 
with sources with very large outer lobes. In this scenario, 



the triple source B 0925+420 ( |Brocksopp et al.| |2007| ) is an 
example of an object where the jet production is interrupted 
again before the bow shocks in front of the new jets have 
reached the ends of the very large outer lobes. If such jet in- 



objects, such as the recently discovered 4C 02.27 ( Jamrozy, 
|Saikia, fe Konar||2009| ). 
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